Color-spin locking in a self-consistent Dyson-Schwinger approach 
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We investigate the color-spin locked (CSL) phase of spin-one color-superconducting quark matter 
using a truncated Dyson-Schwinger equation for the quark propagator in Landau gauge. Starting 
from the most general parity conserving ansatz allowed by the CSL symmetry, the Dyson-Schwinger 
equation is solved self-consistently and dispersion relations are discussed. We find that chiral symme- 
try is spontaneously broken due to terms which have previously been neglected. As a consequence, 
the excitation spectrum contains only gapped modes even for massless quarks. Moreover, at mod- 
erate chemical potentials the quasiparticle pairing gaps are several times larger than expected from 
extrapolated weak-coupling results. 



I. INTRODUCTION 



In recent years, the possible existence of various color 
superconducting phases in Quantum Chromo Dynam- 
ics (QCD) at low temperatures and sufficiently high 
densities has received much attention (for reviews see 
0, Q). At asymptotic densities, where the problem can 
be treated within a weak-coupling expansion, there is 
compelling evidence [3j that quark matter is in the color- 
flavor locked (CFL) phase where up, down, and strange 
quarks participate in a scalar (spin 0) diquark conden- 
sate [4{. The situation is much less clear at "moderate" 
densities of about 3 to 10 times nuclear matter density, 
which could be reached in the centers of compact stars. 
These densities could be too small to allow for a sizeable 
fraction of strange quarks and, hence, pairing of strange 
quarks with nonstrange quarks might be disfavored. On 
the other hand, the cross-flavor pairing of up and down 
quarks in a two-flavor superconducting (2SC) phase could 
be unfavorable as well if local neutrality and beta equilib- 
rium are taken into account [f| • Indeed, as demonstrated 
in various model calculations, BCS pairing of quarks with 
unequal flavors could be absent in the relevant regions of 
the phase diagram @, 0] . 

In such environments quarks of different flavors might 
pair independently. In this case, if we restrict ourselves to 
pairing channels which are antisymmetric in color, spin-0 
pairing is Pauli-forbidden and the most promising can- 
didates are Cooper pairs with total spin 1. Such phases 
were discussed first in Refs. @, Q and their properties 
were investigated in great detail in Ref. [10] for weak- 
coupling. It was found that, among the various pairing 
patterns, the most stable spin-1 phase is the so-called 
color-spin locked (CSL) phase, where the three spin pro- 
jections are locked to the three color degrees of freedom 
of the antitriplet. This is quite analogous to the B phase 
of superfluid helium-3, where the spin is locked to the 
orbital angular momentum. 

In Ref. the CSL phase has been analyzed within an 
Nambu-Jona-Lasinio (NJL)-type model. In this model, 
the locking of color and spin was introduced through a 



gap matrix 



4>- 



A 7 -J 



(1) 



where 7 = (71,72,73) are the gamma matrices with 
spacelike indices and / = (A7, — A5, A2) are the antisym- 
metric Gell-Mann matrices, acting in color space. It was 
found that for massive quarks all quasiparticle modes are 
gapped with the smallest gap in the spectrum being pro- 
portional to the constituent quark mass. This is an inter- 
esting result because ungapped quark matter in general 
leads to rapid neutron star cooling via the direct URCA 
process which is difficult to reconcile with the empirical 
data. In fact, in a recent analysis of neutron star cooling 
data, it was concluded that in a possible quark core all 
quasiparticle modes should be gapped with the smallest 
gap of the order of 10-100 keV and decreasing with den- 
sity [3. The NJL- model results of Ref. [U turned out to 
be in rough agreement with these constraints. However, 
the size of the smallest gap was found to be very sensi- 
tive to the model parameters (see also Ref. [HJ], where 
a similar model with a nonlocal interaction was used). 
Moreover, in order to obtain gaps which are decreasing 
with density the authors used a density dependent cou- 
pling which had to be introduced by hand. In this situa- 
tion a more fundamental, QCD-based approach is clearly 
desirable. 

A second problem was raised in Ref. Here the 

color-spin locking was introduced in a similar way as in 
Eq. |T]) but with 7 replaced by a transverse gamma ma- 
trix, 7i(p) = 7 — (7 • p) p, where p — p/\p\ and p is the 
3-momentum dependence of the gap matrix. In this case 
it turns out that there is always an ungapped mode, even 
for massive quarks. Of course, the question whether or 
not there are ungapped modes in the spectrum should 
not depend on an arbitrary ansatz for the gap matrix. 
Instead, the gap structure should be the result of a self- 
consistent dynamical calculation searching for the most 
favored ground state. Again, this problem cannot be 
solved within NJL- type models, because the point inter- 
action does not allow for momentum dependent struc- 
tures, like xl. 

Recently, there has been considerable progress in 
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studying QCD at nonzero chemical potential within the 
framework of Dyson-Schwinger equations (DSEs) in Lan- 
dau gauge [H, E3] • Dyson-Schwinger equations pro- 
vide a set of coupled integral equations that are the 
quantum equations of motion of the Green functions of 
a quantum field theory. If they are solved completely, 
the whole dynamics of the theory under consideration 
is understood. Of course, in practice, one has to per- 
form truncations. By now there is a lot of experience 
in setting up a truncation scheme for Landau gauge in 
the vacuum, whose results are in good agreement with 
lattice data (for reviews see [II]). In Refs. [H, E3] this 
scheme has been extended to describe spin-0 color super- 
conductors at nonzero chemical potential. This has been 
done in such a way that at very large chemical poten- 
tials the weak-coupling results of Ref. [lj| are correctly 
reproduced. To that end, medium modifications of the 
gluon propagator through particle-hole excitations have 
been implemented. 

Thus being constrained from "two ends" (lattice cal- 
culations in vacuum and weak-coupling results at very 
large fj.) , the framework has then been applied to the "in- 
teresting" regime of moderate chemical potentials. Here 
the authors found considerable deviations from both the 
naively extrapolated weak-coupling pairing gaps [H| and 
the NJL-model behavior of the dynamical quark masses 

Motivated by these results, we have extended the for- 
malism of Refs. [HI, H3] to describe spin-1 color super- 
conductors in the CSL phase. Main focus is thereby the 
question about the smallest excitation gap, raised by the 
(partially contradicting) results of Refs. [H| and [HI- 
To avoid similar problems related to a specific ansatz 
for the gap function, we set up a fully self-consistent 
scheme, starting from the most general parity conserv- 
ing CSL symmetric ansatz for the normal and anomalous 
quark self-energies. Since in Ref. [H[ (constituent) quark 
masses were found to be important, we also include chi- 
ral symmetry breaking terms. As we will see, these terms 
are in fact crucial and, to our surprise, even survive in 
the chiral limit. 

This work is organized in the following way: In sec- 
tion HD we give a brief presentation of the theoretical 
framework for our self-consistent solution of the Dyson- 
Schwinger equation. The numerical results are presented 
in section Hill In section ITVl we present a physical inter- 
pretation in terms of dispersion relations. We summarize 
what we accomplished and conclude in section fVl 



II. THEORETICAL SETUP 

A. The truncated quark DSE at finite chemical 
potential 

We follow the scheme and notations presented in [u| 
to get a truncated, but closed DSE for the quark propa- 
gator at a finite chemical potential within the Nambu- 



Gor'kov formalism. The renormalized quark DSE in 
Landau gauge with appropriate quark- wave function and 
quark-gluon vertex renormalization constants, Z 2 and 
Zip, respectively, is then given by 



s-\p) 



3»«S (?) 



Z 1F T,(p). 



(2) 



Sq denotes the bare quark propagator, 5 the full quark 
propagator and E the quark self-energy. These objects 
are 2x2 matrices in Nambu-Gor'kov space, consisting of 
normal and anomalous components: 



Mp) 

S(p) = 
E(p) = 



= f s+( P ) 

I s -(p) 

(S+(p) T~(p) 

{t+(p) s-(p) 

E+(p) §"(p) 

$+(p) E"(p) 



(3) 



Inserting this into Eq. ([2]) and formally solving for the 
components of the full propagator, we arrive at 



T =-Z ir (Z 2 S ¥ 1 + Z 1F Y^) <P ± S ± , (4) 
Z\f E d 



s ± 1 = z 2 s t 1 1 7 - v± 



Z\ F ^(z 2 S^ 1 + Z 1f TF\ V 



(5) 



Each of these components are, in principle, matrices in 
flavor, color, and Dirac space. However, since we only 
consider pairing of quarks with equal flavor, the flavor 
structure is trivial, and we can consider each flavor sep- 
arately. 

The bare quark propagator is given by 



s o(p) = {-i^plA ~ W ■ 7 + m o) 1 



(6) 



and Sq(p) = —CSq(—p) t C. Here we have introduced 
the notation uj p = p± + ipt. Note that we have dropped 
possible flavor indices for the quark mass too and the 
chemical potential /i, which are in general flavor depen- 
dent. On the other hand, we exclude "color chemical 
potentials", i.e., constant gluonic background fields, be- 
cause it can be shown by the equations of motion that 
these vanish in the CSL phase. 

Defining the gluon momentum by k = p—q, the normal 
self-energy E + and the gap functions <f> + in our trunca- 
tion are obtained via 



S+b)= ^^/(^I^Aa^+^Aaa^fc 2 ) 



P T 



k 2 + G{k) k 2 + F{k) 



, (7) 



<S>+(p) 



P T 



pL 
(If 



k 2 + G{k) k 2 + F{k) 



(8) 
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Furthermore, we have £~(p) = -C£+(-p) C and 
^ _ (p) = 74 $ + (p) 74: which guarantees the action to 
be real Together with Eqs. (gj) and ©, these ex- 

pressions form a set of integral equations which have to 
be solved self-consistently. 

P^ v and P^ v are the usual projectors in the subspace 
orthogonal to the gluon momentum k^. P^ v is transverse 

to the 3-momentum k, whereas P^ v is longitudinal. The 
functions G(k) and F(k) encode Landau damping and 
Debye screening of the gluons. In our analysis this is not 
included self-consistently, but for simplicity F and G are 
computed using bare quark propagators. The result has 
similarities with the hard-dense-loop approximation at 
asymptotically large densities. (For details, see Ref. fl5j[-) 

The function a s (k 2 ), being referred to as the strong 
running coupling, is the only input in our truncation 
scheme. Analogous to previous investigations of the 2SC 
and CFL phase [TEI, [T3], we employ two different cou- 
plings. One of theses running couplings, named ai(fc 2 ), 
has been determined in DSE studies of the Yang-Mills 
sector 21]. It will serve as a lower bound to the cou- 
pling strength, since it underestimates chiral symmetry 
breaking significantly in our vertex approximation. The 
other running coupling considered is chosen to reproduce 
lattice QCD results for the quark and gluon propagators 
in the vacuum for Landau gauge [HJ and will be labeled 
<^u(k 2 ). In order to vary the renormalization point for 
this coupling, we exploit the multiplicative renormaliz- 
ability of the quark DSE. 

As stated above, taking into account the medium po- 
larization, Debye screening and Landau damping are in- 
cluded. Both, screening and damping of the interaction, 
increase with the interaction strength. As a consequence 
the generated gap functions are relatively insensitive to 
variations of a s (k 2 ). Due to this fortunate instance, our 
approach has quite a predictive power despite the uncer- 
tainties in the effective low-energy quark interaction. 

The renormalization constants are determined in the 
chirally broken vacuum. Inserting the above expressions 
for the normal and anomalous self-energy into the DSE, 
Zip cancels out. We determine the quark wave func- 
tion renormalization constant Zi and the renormalization 
constant Z m , relating the unrenormalized quark mass 
mo (A 2 ) at a cutoff A to the renormalized mass m{v) via 

m (A 2 ) = Z m (v 2 ,A 2 )m(v), (9) 

by requiring 

S Z(P)\ P 2=V* = -ij> + rn{v) (10) 

at a renormalization scale v. The mass dependence of the 
quark wave function renormalization constant -Z 2 turns 
out to be negligible in the considered mass ranges. There- 
fore Zi is determined in the chiral limit. 



B. General CSL symmetric ansatz 

As pointed out before, the truncated Dyson-Schwinger 
equation presented above corresponds to a coupled set of 
integral equations. In order to find a self-consistent so- 
lution, we parameterize the normal and anomalous self- 
energies in terms of scalar functions (so-called dressing 
functions) and determine normal and anomalous propa- 
gators via Eqs. (j4j and j5|). These are then inserted into 
Eqs. and (JSJ to recalculate the self-energies, and the 
whole procedure is iterated until convergence is achieved. 

For this strategy we need to ensure that the basis for 
the parameterization of the self-energies and propagators 
forms a closed set, i.e., that no further structures are gen- 
erated during the iteration. To that end, we must allow 
for the most general parameterization which is consistent 
with all symmetries that are left unbroken in the consid- 
ered phase, i.e., in our case, the CSL phase. 

In the Introduction we loosely characterized the CSL 
phase as spin-1 pairing of quarks with equal flavor, where 
the three spin projections are locked to the three color de- 
grees of freedom. We should now be more precise: In the 
CSL phase, the total (spin + orbital) angular momentum 
of the pair is locked to an SO c (3) subgroup of the SU C (3) 
color symmetry. This means, the CSL ground state re- 
mains invariant under transformations of the quark fields, 
generated by the sum Q 

G = J + I, (11) 

where J = S + L, S = ^75747 is the spin opera- 
tor, L = f x pis the orbital angular momentum, and 
I = (A7, — A5, A2) is the generator of the above-mentioned 
50 c (3) subgroup in color space. One can easily verify 
that the generators G satisfy the commutation relations 

[Gi,Gj] = itijkGk, (12) 

corresponding to a residual SUj+ c (2) symmetry. 

For a CSL invariant ground state, all Green functions 
are also invariant under CSL transformations. Therefore, 
we require for the normal and anomalous propagators 

S + {p,p 4 ) = C/S ,+ (p,p 4 )74C/ t 74, 

T+(p, Pi ) = UT+fapfiC-WC, (13) 

where U S SUj+ c (2). In addition, we restrict ourselves 
to even parity pairing, i.e., we require 

S + (P,P4) = 74<S' + (-p,P4)74, 

T+(p,p 4 ) = -7 4 r+(-p,p 4 )74. (14) 

We will now present a self-consistent parameterization 
under these conditions. Since only quarks with equal 
flavor are pairing in the CSL phase, the different flavors 
decouple completely in our truncation. (In fact, the num- 
ber of flavors only appears in the medium polarization, 
i.e., as a prefactor in the screening functions F(k) and 
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G{k).) As a consequence, we only need to parameterize 
the color-Dirac structure. Nevertheless, the number of 
terms which are consistent with our requirements - and 
therefore have to be considered - is still very large: We 
find that normal and anomalous self-energies as well as 
normal and anomalous propagators must be parameter- 
ized by 20 independent functions each. 

The result can be written in the following form, 

5 

S + (p)= Z 2 J2(- i P-7 S A,i ( p ) ~ * ^P 74 S C,i (p) 
i=l 

+ S bAp) -iliP-lS+^p) ) Pi, 

5 

T+ (p) = Z a £ (i 74 £ • (p) + 74 T+. (p) 

i=l 

+ T+» +ip^T+ ti (p) ) M h 

Z 5 / 

+ E B,iO) -iliP-l^D^P) ) P U 

z 2 5 ( 

Zif i=i v 

(15) 

where Pi and Mj are matrices in color-Dirac space defined 
by 

(Pi)ij = Sij 

(-Fb)ij = -75 (p- 4) 

(Pi)y = 7475 (7-4) 

(-Ri)ij = PiPj 

(POij = niliPj+Pilj) (16) 

and .Mj = 75 j^Pi- The indices i,j are indices in color 
space while the Dirac structure is implicit. As before, 
p = p/\p\. 

This ansatz can be derived by considering the largest 
linear independent basis in color-Dirac space under the 
SUj+c{2) symmetry. It can furthermore be deduced by 
coupling two quark fields in three internal SU (2) rep- 
resentations to an overall singlet. Our notation of fac- 
torizing the right hand sides of Eq. (|15p into four pure 
Dirac terms and the color-Dirac matrices Pi or Mi is 
mainly motivated by the structure of the propagators 
and self-energies in a normal-conducting medium, which 
are unit matrices in color space (= Pi) and are custom- 
arily parameterized by the same Dirac terms. Also note 
that the functions Sp t (p), T^^p) E^(p), and 4>fi(p) 
are renormalization point independent. Moreover, those 
terms with F = A or F = C are chirally symmetric, 



whereas the terms with F = B or F = D break chiral 
symmetry. We want to emphasize, that the tensor struc- 
ture is completely determined by the residual symmetry 
and therefore independent of our truncation. 

The defining equation for the anomalous propagator, 
Eq. (H]), is homogeneous. Therefore we are free to choose 
a global phase for the gap function. In our calculations, 
we chose this phase, such that 

<^>4>d - ^i(- P 4,\p\r (17) 

and therefore, all scalar gap functions are real for p^ = 0. 

Having discussed the most general ansatz for the CSL 
phase under even-parity pairing, we can map out previ- 
ous non-self-consistent ansatze by their specific parame- 
terization of the gap functions as given in Eq. (fTS"]) . For 
the ansatz used in Ref. [Ill ] the mapping is particularly 
simple. It is given by 

0+3 = A (18) 

and all other scalar gap functions vanish. 

A more elaborate ansatz has been chosen in Refs. 0, 
H, Ell • l n addition, there was also a distinction between 
different cases of the CSL phase. One of these cases, the 
so-called "transverse case" maps into our notation by 

#A,2 = 0O,3= 2^ + + ^~)' 

^,3 = ^,2 = |(0 + -0"), (19) 

with (j)^ being the (anti-)quasiparticle gap functions in 
the notation of Refs. @, [TTJ . In weak-coupling calcula- 
tions <j)~~ is set to zero, since it does not contribute to 0+ 
even to subsubleading order Q. 

A second ansatz discussed in Refs. @, [TTJ] is called the 
"mixed case" and breaks chiral symmetry. For this case 
the mapping is 

^,3=4,2 = -<^,2 = !W> + -<n- (20) 

In the weak-coupling limit, where again <fi~ is set to zero, 
this ansatz is energetically disfavored compared to the 
ansatz given by Eq. (fT9|) . 



III. RESULTS 

A. The scalar gap functions as functions of the 
chemical potential 

We begin with the numerical results for the gap func- 
tions (f)^ i at the Fermi surface, i.e., at P4 — and 
\p\ = pp. The Fermi momentum pp corresponds to the 



-1 



0,3 0,4 0,5 



0,6 0,7 

H [GeV] 



0,8 



0,9 



-1 



0,3 0,4 0,5 



0,6 0,7 

li [GeV] 



0,9 




0,3 0,4 0,5 0,6 0,7 0,8 0,9 1 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1 

1^ [GeV] 1^ [GeV] 



Figure 1: The values of the gap functions 0j 2 (upper panel) 
and (lower panel), F = A, B , C, D, at the Fermi sur- 
face for quarks with a renormalized mass of 18.4 MeV at the 
renormalization point v = 2 GeV, obtained with the strong 
running coupling ai(k 2 ). 



momentum where the gap functions have their maxima. 1 
We restrict our discussion to the gap functions </>J 2 and 

4>^ 3 , which correspond to the attractive channels. The 
other gap functions are induced by self-consistency, but 
their size turns out to be negligible. 

Our results as functions of the chemical potential fi 
are shown in Fig. [TJ which has been obtained with the 
(weaker) strong running coupling a\(k 2 ), and in Fig. [2] 
which corresponds to a\\{k 2 ). As one might expect, there 
is a general tendency that the chiral symmetry conserving 



1 In general, the Fermi momenta can be defined via Luttingcr's 
theorem, see Ref. |l5j . In the present case, however, the much 
simpler approximation as being located at the maxima of the 
gap functions works extremely well. As the gaps are relatively 
small, the anomalous propagator is strongly peaked at the Fermi 
surface and, consequently, also the gap function. The error due 
to this approximation is small compared with the shift of the 
Fermi momenta due to the normal self-energy. 



Figure 2: The same as Fig. [T] but for the running coupling 
cvn(A; 2 ) with a renormalized mass of 18 MeV at the renormal- 
ization point v — 2 GeV. 

gap functions, (f>\ i and 4>c i: are more important than 
the the symmetry breaking ones, <i>\ { and (I>d i- One 
important exception is the function (f>~^ 2 which can be 
as big as or even bigger than the symmetry conserving 
functions. On the other hand, cf)\ 2 is relatively small. 

If we compare this with previous works, we see that 
the NJL ansatz, Eq. (|18p . is a rather poor approximation. 
But also the "transverse case" , Eq. (|19l) , misses the effect 
of 2 . Only in the "mixed case" , Eq. ([20]) all important 
gap functions are taken into account, although not as 
independent quantities. In fact, it is interesting to note 
that, within numerical precision, we find that 0<t. 2 = 
(f)\ 3 , in agreement with Eqs. (fTTJj) and PD|) . 2 The other 
identities in these equations, however, are not found to 
hold. 



2 This seems to indicate that this identity is protected by an un- 
derlying symmetry. However, until now, we have not been able 
to identify it. 
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The scalar gap functions are of the order of a few 
MeV. This is about two orders of magnitude larger than 
one would have expected from extrapolating the weak- 
coupling results of the CSL phase to small chemical po- 
tentials. This remarkable result can be traced back to two 
features: First, already for spin-0 phases, our approach 
leads to gap functions which are one order of magnitude 
larger than the extrapolated weak-coupling results 
And second, the suppression of the of the spin-1 gaps 
relative to the spin-0 gap is also about one order of mag- 
nitude less than in the weak-coupling limit. E.g., for the 
"transverse case", the latter is given by a factor of e~ 9 / 2 , 
which is the product of two equal suppression factors of 
e -9/4 ^ o.l related to longitudinal and transverse glu- 
ons, respectively [13]. However, when we switch off the 
longitudinal gluons by hand, we find that their contribu- 
tion is almost negligible in our approach and, therefore, 
the suppression relative to the spin-0 result is only about 
one order of magnitude. 

Comparing Fig.[T]with Fig. [2j we find that the different 
scalar gap functions respond differently to the change of 
the strong running coupling. Typically, the results differ 
by a factor of 2, but the effect is somewhat smaller for 
^ 2 ,^ 3 and 4>~q 3 , whereas for <f>^ 2 it is about a factor 
of 2.5 at \x = 1 GeV and even a factor of 4 at fi = 0.5 GeV. 
This is similar to the quark massfunction in the vacuum, 
which differs by a factor of 2-3 between the two couplings 

m. 



B. The influence of the current quark mass 
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Figure 3: The values of the gap functions 
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'c-2 (upper panel) 

Since quark masses played an important role in the and 4>~b-2 (lower panel) at the Fermi surface for various quark 



NJL- model study of Ref. [llj, we have also investigated masses using the strong running coupling ai(k 2 ). 



the effect of quark masses on the scalar gap functions. In 
Figs. [3] and [4] the values of c^jfc 2 and t/>^ 2 at the Fermi 
surface are shown as functions of the quark chemical po- 
tential fi for various quark masses. Obviously, the gap 
function c/><t 2 is almost independent of the quark mass, 
whereas there is a moderate mass dependence of 2 - 
Indeed, recalling that <\>\ i and c6^ ; break chiral sym- 
metry, it was our expectation that these gap functions 
are induced by nonvanishing quark masses and therefore 
the greater sensitivity of cA^ 2 to the quark mass appears 

rather natural. To our surprise, however, <j)~^ 2 remains fi- 
nite even for vanishing quark masses. This means, chiral 
symmetry is spontaneously broken by the diquark con- 
densates. As we will see in the next section, this has 
important consequences for the quasiparticle excitation 
spectrum. 

The results of Figs. [3] and 2] suggest that the sensitiv- 
ity of (fig 2 to the quark mass is more pronounced for the 
weaker coupling a/(fc 2 ) than for the stronger coupling 
a\i(k 2 ). In order to analyze this behavior more system- 
atically, the value of 2 at the Fermi surface is shown 
in the upper panel of Fig. [5] as a function of the quark 
mass for fixed /i = 400 MeV and both running couplings. 



For the stronger coupling au(k 2 ), we find that 2 is 
almost independent of the quark mass for m q < 1 2 MeV 
and then becomes more sensitive to the mass (dashed 
line). This is quite plausible: At low quark masses, the 
value of 2 is dominantly determined by the sponta- 
neous symmetry breaking, whereas at higher masses the 
explicit breaking plays a bigger role. 

For the weaker coupling a\{k 2 ), on the other hand, 
the spontaneous symmetry breaking part is much smaller 
and we find a sizeable mass dependence right from the 
beginning (solid line). In fact, in this case the mass de- 
pendence becomes weaker at larger values of m q . This 
could be due to the fact that at fixed chemical potential 
the density decreases with the quark mass. Note that 
m q corresponds to the quark mass at the renormalization 
point v — 2 GeV and the effective masses at the scale of 
the chemical potential could be considerably larger. 

The effect of increasing the chemical potential is very 
similar to the effect of increasing the coupling strength. 
This can be seen in the lower panel of Fig. [5] where the 
value of <\>\ 2 & f f ne Fermi surface is shown as a function 
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Figure 4: The same as Fig. [3] but using the running coupling 
an(fc 2 ) and slightly different masses. 



of the quark mass for the coupling a\(k 2 ) and two dif- 
ferent chemical potentials. The dashed line corresponds 
to [i = 400 MeV and is identical to the dashed line in 
the upper panel. Increasing the chemical potential to 
800 MeV (solid line) we find again that the spontaneous 
symmetry breaking becomes much stronger and, hence, 
ci^ 2 becomes less sensitive to the quark mass. 



IV. DISPERSION RELATIONS 

A. Discussion 

Since the CSL phase might be relevant for compact star 
phenomenology, and in particular the question whether 
there is an ungapped mode in the excitation spectrum, 
it will be interesting to see how the dispersion relations 
behave for a self-consistent solution. As mentioned in the 
Introduction, the dispersion relations for massive quarks 
in the CSL phase have been investigated in [ll[ and [li[ 
with contradicting results: The authors of Ref. [l4| used 
a DSE approach in the weak-coupling limit and found 



Figure 5: The values of the gap function 



at the Fermi 



surface as a function of the renormalized quark mass at the 
renormalization point v = 2 GeV. Upper panel: Results for 
fi = 400 MeV and different running couplings. Lower panel: 
Results for the running coupling ai(k 2 ) and different chemical 
potentials. 



that there is a gapless mode. In the NJL-model inves- 



tigation of Ref. [llj, on the other hand, it was found 



that there are only gapped modes in the excitation spec- 
trum and that the smallest effective gap is proportional 
to the constituent quark mass. Since both investigations 
lack self-consistency, we have the chance to gain more in- 
sight into the dispersion relations with our self-consistent 
ansatz. 

To describe the quasiparticle excitation spectrum, in 
principle we have to extract the spectral functions. This 
is a relatively difficult task in our approach because one 
needs to analytically continue the results obtained in Eu- 
clidean space into Minkowski space. As demonstrated in 
Ref. [l(| for the case of spin-0 color superconductors, this 
problem can be solved by applying the maximum entropy 
method (MEM) 

In the following, however, we will pursue a simplified 
analysis, aiming at a more qualitative discussion. To that 
end, we concentrate on the dispersion relations e(|p|) that 
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are given as the solutions of the algebraic equation 

Det[5- 1 (p, P4 )]| jp4=£(pl =0. (21) 

Since our calculations are performed in Euclidean space, 
we can not access the roots of this equation. Therefore we 
apply the following approximations: First, we neglect the 
normal self-energy of the quarks and second, we replace 
the scalar gap functions in < f )± (p) by their values at the 
Fermi surface. 

Despite these simplifications, the resulting dispersion 
relations are extremely long expressions and difficult to 
handle. Therefore we have searched for the solutions nu- 
merically. The results will be presented in the next sub- 
section. 

However, to get a deeper understanding of the under- 
lying mechanisms, it is instructive to perform further ap- 
proximations which allow us to derive analytical expres- 
sions for the dispersion relations. In particular, we wish 
to explore the significance of the gap function <f>g 2 , which 
we found to be non-negligible in the self-consistent solu- 
tions discussed in section IIIII Led by these results, we 
consider a simple ansatz with 

0jl,3 = 0(7,2 = 0C, 0s,2 = <t>B, (22) 

and neglect all other gap structures. The resulting dis- 
persion relations 3 are 

e(\p\) = ±^/(n±\p\f + <f>% (23) 

and 



e(\p\) = ±y ( M ± bl) 2 + \ [4>b ± yffe + 8$. J . (24) 

Hence, the simultaneous presence of the gap <j>g, repre- 
senting opposite chirality pairing, and the gap <f>Q, rep- 
resenting equal chirality pairing, implies that all modes 
are gapped. This is similar to Ref. where the spec- 
trum was found to be fully gapped in the "mixed case" , 
Eq. (12"01) , but not in the "transverse case" , Eq. (|19p . 

B. Numerical results 

After these simplified considerations, we now return to 
the numerical solutions of Eq. (I2ip , using the more com- 
plete, although still approximate, expression for the in- 
verse propagator. Since (pg 2 ^ in the self-consistent so- 
lution of the Dyson-Schwinger equations, we expect that 
all quasiparticle modes will be gapped. 

The resulting dispersion relations at \l = 500 MeV 
for massive quarks with a renormalized mass of 



Note, that in this simplified case pp = fi as we neglect normal 
self-energy contributions and quark masses. 



30 




1 1 1 1 1 1 1 1 1 1 1 

0.48 0.49 0.5 0.51 0.52 

|p|[GeV] 

Figure 6: The quasiparticle dispersion relations for quarks 
with a renormalized mass of 18 MeV at the renormalization 
point v — 2 GeV obtained with the strong running coupling 
a u (k 2 ) at fx = 500 MeV. 
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Figure 7: The smallest quasiparticle energy gaps for massive 
quarks obtained with the running coupling a\{k 2 ) (top) and 
cvn(fe 2 ) (bottom). 
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m q {2 GeV) = 18 MeV are shown in Fig.[6j As expected, 
we observe only gapped modes in the excitation spec- 
trum. Of special interest for compact star cooling phe- 
nomenology is the behavior of the size of the smallest 
effective gap as a function of chemical potential. This is 
shown in Fig. [7] for the effective strong running couplings 
a\{k 2 ) and a\\{k 2 ) and for various quark masses. Since 
we find </>£ 2 to be nonvanishing in the chiral limit, cor- 
responding to a spontaneous symmetry breaking, there 
are only gapped modes in the excitation spectrum, even 
for vanishing quark mass . In addition we observe that 
the dependence of the effective gap on the quark masses 
is weak, particularly for higher values of the chemical 
potential. Although the value of the smallest gap is sen- 
sitive to the coupling, we can conclude, that the gap- 
function is monotonically increasing with the chemical 
potentials and is of the order of a few 100 keV in the 
astrophysically relevant regime. 

V. SUMMARY AND CONCLUSIONS 

We have studied the quark propagator in the color- 
superconducting CSL phase at zero temperature in a 
self-consistent Dyson-Schwinger approach. Our main ob- 
jective has been the clarification of contradicting pre- 
dictions by previous investigations, especially concerning 
the phenomenologically interesting size of the smallest 
excitation gap. To this end, we have considered the most 
general ansatz for the quark propagator allowed by the 
CSL symmetry pattern. For comparison, we have out- 
lined the truncated ansatze of previous investigations and 



how they are parameterized in our basis. In our numeri- 
cal exploration, we found that the CSL phase exhibits a 
color-spin structure that is close to the color-spin struc- 
ture used in weak-coupling calculations [lOj. As a re- 
markable difference, chiral symmetry breaking terms are 
additionally dynamically generated, leading to a spon- 
taneous symmetry breakdown. As a result, in contrast 
to [14|, the size of the smallest excitation gap is non- 
vanishing and triggered by a chiral symmetry breaking 
diquark condensate. Therefore, the dynamical mass gen- 
eration, as elaborated in (TT| , is not necessary for obtain- 
ing a nonvanishing smallest excitation gap. Although the 
value of the smallest gap is sensitive to the coupling, we 
can conclude, that the gapfunction is increasing mono- 
tonically with the chemical potentials and is of the order 
of a few 100 keV in the astrophysically relevant regime. 
This seems to be in conflict with the analysis of Ref. [l2j 
where it was concluded that the smallest gap should de- 
crease with the chemical potential in order to describe the 
neutron star cooling data [12]. However, before drawing 
strong conclusions, a more careful analysis is certainly 
required, where both specific heat and neutrino emissiv- 
ities, are derived within our approach. 
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